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Abstract 

We prove that the ground state momentum distribution of a one-dimensional system of impene- 
trable bosons exhibits a k~ 4 tail for any confining potential. We also derive an expression for easily 
computing the asymptotic occupation numbers and verify our results with an exact numerical 
approach. 
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I. INTRODUCTION 

It is well known, since Girardeau original paper that one- dimensional impenetrable 
point-like bosons can be mapped into a free fermion problem taking the product of the many- 
body wave function of the fermion system with a unit antisymmetric function. Although 
restricted to one dimension, this remarkable observation opened the possibility of treating 
analytically a strongly interacting problem. As can be easily seen local observables (such 
as the probability of finding the particles in a well defined position) remain the same in 
both systems, but non-local quantities of the bosonic problem need to be calculated in a 
non-trivial way. In this letter we focus on the one-body reduced density matrix n(x, y) for 
an iV-body system of bosons at T = 0: 

n(x, y) = N J tp*(x, x 2 , . . . , x N ) ij)(y, x 2 , ■ . ■ , x N ) dx 2 , . . . , dx N , (1) 

where ib is the many-body wave function for the ground state (unit norm). As pointed 
out in [l| it is not an easy task to reveal the properties of this object because the multiple 
integration can not be carried out explicitly. Difficulties remain in quantities such as the 
momentum distribution (MD), which can be found projecting equation ([T|) onto a plane 
wave basis. The low momentum region of the MD was explored by Schultz and Lenard 
to determine whether Bose-Einstein condensation occur in the thermodynamic limit, 
as argued by Girardeau {if. However, at that time, the interest in the model was mainly 
theoretical and two decades pass until further analytical results concerning the asymptotic 
properties of the momentum distribution were published 4|. 

This situation changed ten years after the first successful experiments on Bose-Einstein 
condensation of alkali gases 51. These last provided the framework for achieving effectively 

n □ ft 

one-dimensional (T « 0) bosonic systems [6( in laboratory experiments [3, |8[. After its 
experimental realization, further properties of the momentum distribution were extensively 
studied (mainly numerically) not only for theoretical reasons but because the MD char- 
acteristic features helped to determine whether the system had entered into the strongly 
interacting regime. 

In the following we show an asymptotic expansion for the one-body reduced density matrix 
valid for short-range correlations and we derive an universal power law for the tail of 
the momentum distribution, namely k~ A . This behaviour was found previously only for 



particular potentials 



ll| or constrained to formulations on the lattice 12J 



The goal of this paper is to give a general and self-contained proof of this power 
law valid for any confining potential. In addition, we derive an expression which allows 
an easy computation of the short distance properties of the one-body reduced density matrix. 



II. SHORT DISTANCE CORRELATIONS 



Consider a one-dimensional system o 
potential V(x). As shown by Girardeau 
be constructed as: 



N impenetrable point-like bosons trapped by a 
1] the ground state many-body wave function can 



V>(xi, . . . , x N ) = —f= A(xi, . . . , x N ) det [<j>i{xj)], (2) 

where A(x±, . . . ,xn) is the unit antisymmetric function and {(f>i}f =1 are the N lowest eigen- 
states of the one-body Hamiltonian: 

E l U*) = -^%+V{x)Ux), (3) 

which can be taken as real functions. To simplify our notation we will continue to consider 
the case of the ground state but the same arguments will hold, with minor modifications, 
for any other pure state of the system (see Section [TTTT) . 
Definition (JTJ can be written as: 

n(x, y) = — ^ J sgn(x - x 2 ) . . . sgn(x - x N ) sgn(y - x 2 ) . . . x 
. . . sgn(y - x N ) I ^ Sgn(a) ct( i) (x) ct(2 ) (x 2 ) . . . <f) a (N) (%n) ) x 

Sgn(X) (j>x(i)(y) (t>\(2){x 2 ) ■ ■ ■ <P<t(n)(zn) J dx 2 , ■ ■ ■ , dx N , (4) 

where sgn(x) stands for the unit antisymmetric function and Sgn(a) for the sign of the 
permutation a G Sat. This equation can be rearranged in the following form (assuming, 
with no loss of generality, that y > x): 

n(x,y) = , N _ 1 y Yl ^(i)(^)0A(i)(y) x 
Sgn(a) Sgn(X) f[ (s<r(j),m ~ 2 f ^{j)( z )4>x{j)(z) dzj . (5) 



Equivalently one can write an expression removing the constraint y > x, namely: 

n(x,y) = <f><T(i)( x )<f>m(y) S 9n(o-)Sgn(\) x 



(6) 



where #(x) is the Heaviside step function. This expression is exact an can be used to cal- 
culate n(x,y) numerically when the eigenfunctions are known (or calculable). Note that 
the product of all Kronecker-delta terms corresponds to a free fermion distribution and the 
remaining "bosonic correction" appears with powers of J y . To proceed further we note that 
short distance properties of this expression should be well approximated to first order in 
f , because the projection of .urge scales will be strong* suppressed in the high momen- 
turn regime (observation stated without proof in [13|). This is equivalent to saying that 
occupation numbers on a basis {Xi}i>i ( n °t necessarily an eigenbasis of the Hamiltonian) , 
(hk) = J x*k{ x ) n ( x i y)Xk(y)dxdy, will be well approximated by the first order term of n(x, y) 
if k » N (N fixes a maximum energy scale in the system). As we shall see in the last 
section these facts are correct, the leading correction is of first order in f and the remaining 
terms with more than one integral can be neglected. We will first explore the consequences 
of this argument and leave the details of the proof to the last section. 
As claimed before, performing the sum up to first order in JJ yields an approximate expres- 
sion (n s d) valid for describing the short distance properties of the system (y > x): 

N N 

n sd (x,y) = Y^faWMv) ~ 2 / ( Pi{ x ) ( t ) i{y) ( l ) j{ z ) ( t>j{ z ) dz + 
i=l i,j=l Jx 

N .y 

+2^ / faWfatyWiWfaWdz . (7) 
ij=l Jx 

Note that no permutations are involved in this sum, which makes it easier to handle. As 
an example of the applicability of this expression to study short distance properties let 
us compute the MD for a flat confining potential, which has been extensively studied in 
literature. The approximate expression for the occupation numbers in a box of unit length 
can be evaluated as: (n k ) sd = J xl( x ) n sd(x, y)Xk(y) dxdy, where Xk( x ) = {v / 2sin(7r/ox)} fc >i. 
This expression involves simple combinations of trigonometric functions and can be evaluated 
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explicitly. We have carried out this calculation finding: 

N 



2kH 2 -k 2 f -i 2 f 



7T" 



(8) 



As a first check of our approximation we have taken the thermodynamic limit of equation 
(JSJ), finding (h k ) « ^(ir) 4 - ^his expression agrees with previous calculations after the 
identification of N with the Fermi level kw 11]. 

The convergence of the approximated MD to its exact value is rather fast. We have verified 
this property comparing the results of the approximation with the exact ones (equation (jSJ) 
or ([6])) for a finite number of particles. As a concrete example we computed the exact mo- 
mentum distribution for a system of 5 bosons in the flat one-dimensional box. We exhibit 
these results in Figured] where we plot the exact numerical evaluation (hk) and the approxi- 
mate one (hk) sd - As argued before, when k >> N, the agreement between both expressions 
is excellent. 
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FIG. 1: Numerical evaluation of the exact distribution for 5 bosons in a flat box (hj.). Comparison 
with the approximated short-distance expansion (n^.) S(i (from equation ([8])). Note that both axes 
are dimensionless. 
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Finally we note that equation (|Sj) can be approximated for large k by an expression of 
the form k~ 4 + 0(k~ 6 ). In the next section we show that this is a general characteristic of 
a confined Tonks- Girardeau (TG) gas. 



III. ASYMPTOTIC TAIL 

Let us now show that the momentum distribution decays as k~ 4 regardless of the confining 
potential. For a confined system of N hard-core bosons in a general potential V(x) take 
{4>i}f =1 to be the set of the N lowest energy eigenfunctions. As argued previously the 
occupation numbers for large momentum k can be determined projecting equation ((Tj) for n sd 
onto a plane wave basis. The projection of the first term in that equation can be disregarded 
in almost all situations. In fact, if the potential is a smooth function, say C a piecewise, then 
any eigenfunction is at least C a+2 . This means that the coefficients of the projection will be 
of order l/k 2( * a+2 \ So, for smooth functions V(x) the asymptotic behaviour is determined 
by the second and third terms. From now on we omit the first term in equation ([Tj). Then, 
coming back to: 

(™fc) s <2 = J Xl(x)n sd (x, y)xk{y) dxdy , (9) 
we explicitly remove the condition x < y from the expressions. In order to do this we write: 

n sd (x,y) = N(x,y)[9(y-x)-9(x-y)} , (10) 

where Vx, y 

N(x,y) = -2 J2 / Mx)My)Mz)M*)dz + 

+2^ ! faWfaivWiWhWdz . (11) 
i,j Jx 

To perform the integral ([9]) in a plane wave basis Xk = we change variables to: 

(12) 




Now we have: 



J e ikr Hjr J N(X,r)dX^j [0(r) - 6{-r)]dr . (13) 



The power law for the high momentum regime can be found expanding N(X, r) in Taylor 
series around r = 0. Surprisingly enough a direct calculation shows that the function N(X, r) 
verifies: 

N(X,0) = 

f(*,o) = o 

(14) 

0(*,o) = o 
k 0(x,o)^o 

Which is a consequence of its symmetries. Taking (IT41) into account, one ends with (up to 
leading order in r): 

(n fc ) sd oc J e ifc V 3 [0(r) - 0(-r)]dr = ^ / e^[0(r) - 6{-r)]dr. (15) 
Which is equivalent to (n k ) sd oc &~ 4 . 

Notice that the same argument holds when excited states of the gas are considered. In 
such situation the set of energy eigenfunctions of the one-body problem which are used 
to calculate the many-body wave function will have one element which corresponds to the 
largest eigenenergy. This element fixes the scale ko (maximum typical wave number involved 
in that energy scale) such that equation (JTj) becomes a good approximation in the region 
k ^> k (in the numerical example of the previous section k = N). 

IV. TECHNICAL DETAILS 

In the previous section it was shown that the leading order contribution of the term 
with one integral is of order r 3 . As a consequence of this fact the approximate expression 
(7) would be valid only if we prove that the terms with two and three integrals vanish in 
equation ([6]), which is a non-trivial fact (note that those terms can, in principle, contribute 
to order r 3 ). To see this, we first note that the terms with two integrals have a prefactor 
[9(y — x) — 6{x — y)] 2 = 1, which is non-singular and thus can be neglected as we argued for 
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the fermion term in equation (9). Secondly, for terms with three integrals one gets: 

( t ) ^w( x ) ( l ) Hi)(y)[ e (y- x )- e ( x -y)] s 9 n (^)S9n(x) x 

(Pa(a)(z)(f)x(a)(z) dz^j (J <f> a (b){z)<f>\(b){z) dz\ X 

( I (j)a(c)(z)(j)x(c)(z) dz) f[ S aU)Mj)j (16) 

where the complete term is the sum over all choices of the integers a, b and c. For each fixed 
choice the contribution to the order r 3 is found by replacing each integral to first order in 
r. The resulting prefactor reads: 

Sgn(a)Sgn(X)cf> a{1) (X)4> x{1) (X)x 

CT (a)(X)0 A(a) (X)0 (T(6) (X)0 A(6) (X)0 (T(c) (X)0 A(c) (X) , (17) 

which vanishes when summing over all permutations. This completes the proof. 



V. CONCLUDING REMARKS 

We have shown that the momentum distribution of a TG gas in the high momentum 
limit can be described by a power law of the form k~ A regardless of the confining potential. 
The insight given by equation (J7|) shows that this universal behaviour is a consequence of 
the wave function symmetries together with the cusp condition imposed by the subsidiary 
requirement i/j(xi, . . . , x n ) = if = Xj. 

We have also proven the validity of equation (J7J), which is found to be a suitable alternative 
for an accurate and numerically-efficient description of the high-momentum region of the 
TG gas. 
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